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In this talk I will discuss some of the techniques that have been developed over the 
past 35 years to estimate the energy of charged matter. These techniques have been 
used to solve stability of (fermionic) matter in different contexts, and to control 
the instability of charged bosonic matter. The final goal will be to indicate how 
these techniques with certain improvements can be used to prove Dyson's 1967 
(-H ' conjecture for the energy of a charged Bose gas — the sharp TV 7 / 5 law. 
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1. INTRODUCTION 

It is my aim in this contribution to review the main techniques developed to study the 
problem of stability of matter or rather stability or instability of ordinary matter. By or- 
dinary matter I mean a macroscopic collection of charged particles (nuclei and electrons) 
interacting solely through electrostatic or electromagnetic forces. 

The problem of stability of (ordinary) matter interacting through electrostatic Coulomb 
forces was first formulated by Fisher and Ruelle |18j and proved in the seminal papers of 
Dyson and Lenard |12l although, of course the problem of stability of individual atoms 
goes back to the origin of quantum mechanics. An important assumption needed is that 
either the negatively or the positively charged particles are fermions, i.e., obey the Pauli 
exclusion principle. Without the fermionic assumption there is no stability as proved by 
Dyson in 

The importance of the Pauli exclusion principle for stability had been pointed out in 
the celebrated work of Chandrasekhar |H] on gravitational collapse and stability of white 
dwarfs. It was, however, not until the work of Dyson and Lenard that the importance of 
the exclusion principle for the stability of ordinary matter was emphasized. In jl 1 j Dyson 
makes a very precise conjecture regarding the nature of the instability without the exclusion 
principle. 
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It is my ultimate goal here to discuss the proof of this conjecture, which for the main 
part, is joint work with E.H. Lieb. 

Since the work of Dyson and Lenard there has been a lot of activity in the area of stability 
of matter. Most celebrated is the work of Lieb and Thirring [2] giving an elegant proof 
with a bound of the correct order of magnitude. Several variations of the model have been 
studied. Rclativistic effects and magnetic interactions have been included. I shall review 
some of these results below. In recent years the attention has turned to studying the effect 
of quantizing the electromagnetic field, what is often referred to as non-relativistic quantum 
electrodynamics. I will not get into this recent development here. 

2. CHARGED MATTER IN QUANTUM MECHANICS 

The systems we discuss here are all given by iV-particle Hamiltonians (N being a positive 
integer) of the form 

3=1 l<i<j<N 1 31 

Here Zj € K is the charge of particle j, Xj € M 3 is the coordinate of particle j, and Tj is the 
kinetic energy operator for particle j. We shall here consider the kinetic operator to be of 
one of the following four types. 

(1) The standard non-relativistic kinetic energy 

Ti = —A,. 

2m, J 

(2) The relativistic kinetic energy 

If el = yJ-JAj + mjc 1 - m 3 c 2 . 

(3) The magnetic kinetic energy 

2f* = J- f-<V j-^A( Xj 
3 2m j V J c J 

(4) The magnetic Pauli kinetic energy 

yPauli = 1 / ( _ iV _ *LA( Xj )) ■ <T. 

3 2mj V c 
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Above, rrij > refers to the mass of particle j, A : M 3 — ► M 3 is the vector potential for the 
magnetic field, and <jj is the vector of 2 x 2 Pauli matrices for particle j, i.e, er = (ai, 02, 03), 
with 




Ol= , CT2 = , CT3 





The Pauli kinetic energy acts on spinor valued functions (spin 1/2 particles), see the dis- 
cussion of the relevant Hilbert spaces below. 
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The term U in Q represents the magnetic field energy 



U= ^- [ |V x A{x)\ 2 dx 



We are using units in which Planck's constant h — 1 and the fundamental unit of charge 
e = 1. In these units the physical value of the speed of light c is approximately 137 or more 
precisely the reciprocal of the fine structure constant a. 

Many objections may of course be raised about the Hamiltonian H^. Even with the 
relativistic kinetic energy operator the Hamiltonian is not truly relativistically invariant. 
For spin 1/2 fermions the appropriate kinetic energy operator would rather be the Dirac 
operator. Also the magnetic field is treated classically as opposed to as a quantized field. 
Unfortunately, we do not today know how to describe a mathematically consistent fully 
relativistic model in which to even formulate the problem of stability. As is well known 
the Dirac operator is not bounded from below and introducing quantized fields requires a 
renormalization scheme. These questions are being actively studied at present, but we shall, 
as already mentioned, not discuss them here (see One approach to the Dirac 

operator would be to restrict to positive energy solutions for the free Dirac operator. For 
the non-magnetic case this would lead to a problem essentially identical to the problem 
with the relativistic kinetic energy T Rcl . In the magnetic field case this is somewhat more 
difficult (see [27]). 

We now turn to the important question of which Hilbert Space the operator Hn acts 
on. The operator will be unbounded, but we shall always consider the Friedrichs' extension 
of the restriction to C°° functions with compact support. Since the problem of stability is 
a question of lower bounds, the Friedrichs' extension is the correct setting. 

We begin by discussing the Hilbert spaces for just one particle, the one-particle space. 
We consider, in general the one particle spaces to be the square integrable functions corre- 
sponding to particles with q internal states, i.e., Tif = L 2 (R 3 ;C 9 ), where the non-negative 
integer q may be different for the different particles. In terms of spin, this corresponds to 
particles of spin (q — l)/2. 

The many particle space of interest is then 



where K is a positive integer. In order for the Hamiltonian Hn to act on the above Hilbert 
space we must require that all the masses and charges of the particles j = 1 , . . . , N — K are 
the same. Put differently, these particles are identical fermions. For simplicity we assume 
that irij = 1 and Zj = —1, for all j = 1, . . . , N — K, (Except for the sign of charge this 
simply amounts to a choice of units). Thus the fermions are negatively charged. We assume 
the remaining particles to be positively charged, i.e., Zj > for j = N — K + 1, . . . , N. 
Moreover, all the fermions must be described by the same kinetic energy operator. The 
positively charged particles could, in principle, have any kinetic energy operators and even 
simply vanishing kinetic energy (equivalent to infinite mass). Of course, the Pauli kinetic 
energy operator requires the spin of the particle to be 1/2, i.e., q = 2. 

Since fermions have half-integer spin one would maybe like to restrict the negatively 
charged fermions accordingly, i.e., the corresponding q to be even. For the discussion here 
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this restriction, however, plays no role. 

The reader may ask why one does not consider the situation when all or some of the 
positively charged particles are fermions or bosons. In fact, stability in this case is a simple 
consequence of the situation discussed above, since considering fermions or bosons simply 
means restricting to certain subspaces. 

To study what happens when we ignore Fermi statistics all together we shall consider 
all particles to be described by the standard non-relativistic kinetic energy T and all to 
have mass mj — 1. Moreover, we assume that all the charges Zj are either +1 or — 1. We 
indeed consider the charge of each particle to be a variable. Put differently, we consider the 
iV-particle Hilbert space 



n 



N 



lL 2 (M 3 x {1,-1}), 



(3) 



where the set {1,-1} refers to the charge variable. 

We define the ground state energy as the bottom of the spectrum of the Hamiltonian. 
We will not here address the question of whether this is and eigenvalue or not, i.e., whether 
there really is a ground state. When a magnetic field is present we consider the vector 
potential as a dynamic variable and we also minimize over it. Thus 



E 



Physical 
N ' 



inf inf spea^phys 



a !(*,*) 



physical 
N 



and 



E N = inf spec Hn H 



N 



^c fl M n^\ {0} 



\{0} (4) 



(5) 



Here C§° refers to smooth functions of compact support. In we minimize over all smooth 
vector potentials for which |V x A\ E L 2 (R 3 ). Of course the minimization over the vector 
potential is only relevant if the kinetic energy really depends on A otherwise the minimum 
simply occurs for A = 0. 

The energy Em on the space where we consider the charges as variables is precisely the 
same as we would get if we calculate the energy for any fixed choice of charges and afterward 
minimize over this choice. Note, in particular, that the charge variable commutes with the 
Hamiltonian H^. The reason for including the charges as variables is that the Hamiltonian 
is then fully symmetric in all N particles and not just in the positively or negatively charged 
particles separately. 

Strictly speaking the energies in (@J and (jSJ are only defined as the bottom of the 
spectrum if the rightmost expressions are finite (i.e., not — oo). It is in this case that we can 
define the operators as Friedrichs' extensions. The property that the ground state energies 
are bounded below is often referred to as stability of the first kind. It holds except in the 
cases when the fermions are described by the kinetic energy operators T Rel or j ,Pauh . I n 
these two cases stability of the first kind requires that max{zj}/c — maxj{zj}a is small 
enough in the case of T Rcl or that qmaxj{zj} /c 2 = qma.Xj{zj}a 2 is small enough in the 
case of T Pauli . 
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3. STABILITY AND INSTABILITY OF MATTER 

Stability of matter is a stronger statement than stability of the first kind. It means that 
there exists a constant C £ R such that 

^Physical > _ CN (6) 

for all N, i.e., that the total binding energy per particle is bounded. 
Theorem 3.1. (Stability of Matter) 

On the space 7i^ r hyslcal stability of matter 0) holds with a constant C that depends only on 
q and maxjzj} if anyone of the following situations hold. 

• The fermions are described by the standard non-relativistic kinetic energy T or by 
the magnetic kinetic energy T Mag . 

• The fermions are described by the relativistic kinetic energy T Rel and qa(= qc^ 1 ) is 
small enough and maxj{zj}a < 2/ir. 

• The fermions are described by the Magnetic Pauli kinetic energy T Pauh and qa{= 
qc^ 1 ) and qmax.j{zj}a 2 are small enough. 

The case of the standard non-relativistic kinetic energy is the situation first settled by 
Dyson and Lenard j!3j and later by Lieb and Thirring 3j with a constant of the correct 
order of magnitude. There was also a proof by Federbush |14| . 

Stability for the magnetic kinetic energy is an immediate consequence of the non- 
magnetic case and the diamagnetic inequality. In fact, the ground state energy is achieved 
without a magnetic field. 

Stability for the relativistic kinetic energy was first solved by Conlon j^] and then im- 
proved by Fefferman and de La Llave |17| . The version formulated here (which is sharp 
with respect to the bound 2/ir) is due to Lieb and Yau [22] • The case of one electron and 
one nucleus had been studied previously by Herbst [22] and Weder PHj and the case of one 
electron and several nuclei by Lieb and Daubechies 

A proof of stability for the magnetic Pauli kinetic energy was first published by Lieb, 
Loss, and Solovej |2fi| . but had been previously announced, (although with weaker bounds) 
by Fefferman (only later published in |15|). The result as formulated here is optimal in the 
sense that if either qa or q maxj{,Zj}a 2 are large then stability does not hold. This was 
realized in a series of papers |2()l 1241 13"3"| . 

In the case without Fermi statistics Dyson JT] proved that there is no stability and he, 
in fact, conjectured the following result. 

Theorem 3.2. (The asymptotic energy of a charged Bose gas) 

For the energy En defined in 0) we have the asymptotics 



<j>o, <t? = iy <7) 



where 



J= - / l + x A -x 2 (x A + 2y' z dx 
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The reader may wonder why the theorem refers to a charged Bose gas when, in fact, no 
statistics was enforced in the definition of the Hilbert space Hn- The reason is that the 
ground state energy on Hn is the same as the ground state energy one would get if re- 
stricted to the fully symmetric subspace (i.e., the bosonic subspace). This is a fairly simple 
consequence of the facts that the Hamiltonian is fully symmetric in all variables and that 
the expected energy of any trial state does not increase if we replace it by its absolute value. 

It is a simple consequence of the classical Sobolev inequality (see below) that the right 
side of J7J is finite. Dyson proved in JT] that En < —CN 7 / 5 , but not with the correct 
constant. His method and conjecture was inspired by the Bogolubov approximation, which 
had been previously used by Foldy ^5] to calculate the energy asymptotics for the high 
density one-component Bose plasma (bosonic jellium). The Bogolubov approximation is 
usually applied to bosonic systems, it is therefore important that we can think of our system 
as such. 

In their original work on stability of matter |12j Dyson and Lenard obtained as a corollary 
a lower bound on the Bose energy En > ~C'N 5 ^ 3 (see also Brydges and Federbush 
The correct exponent 7/5 was however first proved by Conlon, Lieb and Yau in [?], where 
one finds the lower bound En > —C"N 7 / 5 , but again not with the correct constant. 

Dyson's conjecture was finally proved in two papers establishing respectively an upper 
and a lower bound. An asymptotic lower bound of the form (|7J) was proved by Lieb and 
Solovej |3()| . The corresponding asymptotic upper bound will appear in Solovej |34| . 

The rigorous calculation of Foldy's high density asymptotics for the energy of the one- 
component plasma can be found in Lieb and Solovej (lower bound) and Solovej jHU 
(upper bound). 

The main goal here is to review the techniques used to prove Dyson's conjecture and 
to sketch the main steps in the proof. Many of the techniques used were developed to 
study stability of matter and I will therefore use the opportunity to briefly review these 
applications as well. 

The main strategy in proving the lower bounds in Theorems 13.11 and 13 . 21 is to estimate 
the Hamiltonian below by a simplified Hamiltonian for which one can get a fairly explicit 
lower bound. For Theorem 13 . II the simplified Hamiltonian is a non-interacting mean field 
type Hamiltonian. For Theorem 13.21 the simplified Hamiltonian is a non-particle conserving 
Hamiltonian, which in the language of second quantization is quadratic in creation and 
annihilation operators (a quadratic Hamiltonian). 

We first discuss how to treat the simplified Hamiltonians. 



4. ONE-PARTICLE OPERATORS AND QUADRATIC 
HAMILTONIANS 

The most basic estimate on the ground state energy of a one-particle Hamiltonian, given in 
the next theorem, is a simple consequence of the Sobolev inequality J \ Vip\ 2 > C(J j^l 6 ) 1 / 3 , 
C > for functions ip on R 3 . 

Theorem 4.1. (Sobolev estimate) 

Let V be a locally integrable function on R 3 , such that the Schrddinger operator — A + V 
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may be defined as a quadratic form on C°° functions with compact support. Then 

-\a + V>-C s J\V{x)fl 2 dx 
where Cs > and |t|_ = max{— 1,0}. 

This theorem of course immediately implies that the many-body operator 

N 



has the lower bound —CgN J \V(x)\ 5 / 2 dx. On the fermionic subspace we, however, have 
the much stronger Lieb-Thirring inequality |31| . 

Theorem 4.2. (Lieb-Thirring inequality) 

On the fermionic space /\ M TL\ (M a positive integer) we have the operator inequality 

M 



J2~\ A * + V ( x i) > -Clt? J \V(x)\ 5 J 2 dx, 



where Clt > V is a locally integrable function on R 3 . 

For the relativistic operator T Rel and the magnetic Pauli operator T Pauh one has similar 
inequalities on the space f\ M H\ (q = 2 for T Pauli ). For T Rcl (with m = 1) Daubechies 
proved 

M . 

?f 01 + V(xi) > -C B qc- 3 / |V(a;)|id!B - Mc 2 . (8) 
i=i •* 

For T Pauh (with m — 1 and z — 1) the inequality 

M 



<TTf^ + V{ Xi ) > -C$ s J \V(x)fJ 2 dx 



3/4 / r \ 1/4 



-C^ s c' 3/2 / |V x Afdx / \V(x)\idx ) , (9) 



can be found in |26j . 

When studying bosonic systems one may, as explained above, use the Sobolev estimate 
to get a lower bound on the energy. This will in general not give the best dependence on the 
number of particles. E.g., for the charged gas this would lead to a lower bound — C7V 5 / 3 as 
in [121 ITT] . To get the sharp dependence —CN 7 / 5 a more precise treatment of the interplay 
between the kinetic energy and the Coulomb potential is needed. 

This brings us to Bogolubov's method for calculating the energy of a Bose gas. In the 
Bogolubov approximation one assumes that most particles form a Bose condensate in a 
momentum zero state. The main contribution to the ground state energy will come from 
the correlation between two non-condensed particles with opposite momenta. This effect 
is most easily explained using the formalism of creation and annihilation operators. The 
following theorem gives a rigorous formulation of Bogolubov's method in the simple case 
used by Foldy in |19j . 



February 7, 2008 14:39 WSPC/Trim Size: lOin x 7in for Proceedings solovej-icmp-june3 



8 

Theorem 4.3. (Bogolubov's method) 

Assume that b± ± are four (unbounded) operators defined on a common dense domain on 
a Hilbert space, such that their adjoints are also defined on the same domain. Assume 
moreover that in the sense of quadratic forms on this domain we have the commutator 
identities 

[b* T ,, z ,X.A = - [K>,-X,+] = [K>,+>K,-] =0. forallz,z',T,T' = ±, 

and 

[b TZ , b* J < 1. for all z,t = ±. 
Then for all t,g+,g- > we have (again in the sense of quadratic forms) 

r,z— ±1 z,z'=±l 

> -(* + g+ + g-) + V(t + "+ + 9-) 2 -(g+ + g-) 2 - 

The operator inequality in this theorem follows simply by completing squares (see also |30|1. 
The importance of the lower bound is emphasized by the fact that it is sharp if the operators 
&±,± are truly annihilation operators. In Bogolubov's method the operators are however 
not exactly annihilation operators. Rather one should think of b± tZ as the operator that 
annihilates a particle with charge z and momentum ± (here only the sign of the momentum 
is important) and creates a particle in the condensate (i.e., with momentum 0) with charge 
z. The value t represents the kinetic energy of the particles with momentum ± and g± 
represent the strength of the Coulomb interaction. 



5. REDUCTION TO A SIMPLIFIED HAMILTONIAN 

The reduction to a simplified Hamiltonian require controlling the Coulomb interaction. It 
is often convenient to replace the Coulomb potential by the Yukawa potential 

\x\ 

Replacing the Coulomb potential lo by a Yukawa potential with [i > amounts to 
introducing a long distance cut-off in the potential. It is easy to control this replacement 
since Yq — Y^ has positive Fourier transform (is of positive type) . Hence 

JV 2 

l<i<j<N i=l 

since (Y„- Y M )(0) = fi. 

In order to bound the many-body Hamiltonian Hn below by a one-body Hamiltonian 
one must estimate the two-body potential by a one-body potential. 

Theorem 5.1. (Onsager's correlation estimate) 

Given particles with with charges Z\, . . . , Zn at positions x\, . . . , xjq ■ Let 

Di = min{|a; l — Xj\ \ z.Zj < 0}, 
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i.e. Di is the shortest distance from the particle at Xi to a particle with the opposite charge. 



Above Di depends only on Xi and on the positions of all the particles with the opposite 
charge of the particle at xi. This theorem allows us to separate the original Hamiltonian 
Hn in a one-body Hamiltonian for all the negatively charged particles (with a one-body 
potential depending on the positions of all the positively charged particles) and a one-body 
Hamiltonian for all the positively charged particles (with a one-body potential depending 
on the positions of all the negatively charged particles). 

The proof of this theorem essentially goes back to Onsager (for fi = 0), who addressed 
a stability question for classical charged matter. In [3D] it is being used to introduce a short 
distance cutoff in the potential (i.e., to replace Yq by Y"o — Y M f° r l ar g e A*)- The one-body 
potential is controlled using the Sobolev estimate Theorem 14. II 

If both the negatively and positively charged particles are fermions, we may use Onsager's 
estimate together with the Lieb-Thirring inequality in Theorem 14.21 to prove stability of 
matter. This special case was treated by Dyson and Lenard in their first paper |12| . where 
they also used a version of Onsager's correlation estimate, but did not have the Lieb-Thirring 
inequality at their disposal. 

In order to prove stability of matter as formulated in Theorem 13.11 one may use the 
following stronger correlation estimate due to Baxter £Q . 

Theorem 5.2. (Baxter's correlation estimate) 

Assume as before that all negatively charged particles have charge —1 Then 

N 



The importance here is that the sum on the right is only over negatively charged particles. 
Thus the right side is a one-body potential for the negatively charged particles thinking of 
the positively charged particles as fixed. A sharper (although more complicated to state) 
estimate was given by Lieb and Yau in |.'S2| . The estimate of Lieb and Yau is sharper in the 
sense that the coefficient l + 2z above may be changed to z, but additional bounded errors 
are needed. 

Remark 5.1. (The proof of stability of matter) 

A proof of stability of matter in the standard non-relativistic case may proceed as follows. 
We use Baxter's estimate to arrive at 



Then 





N 



H N > (T i -(l + 2ma^{z. l }){D- 1 -R- 1 ))-N{l + 2msix{z l })R- 1 



i=l,z;<0 



> -CN(l + 2ma,x{ Zl }) 2 , 
where we have used the Lieb-Thirring inequality Theorem 14 . 21 for the potential 

V = (l + 2max{z l })(-Dr 1 +i?- 1 ) 
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and we chose R ~ (1 + 2 max{zi}) _1 . 

The proof of stability of matter by Lieb and ThirringPT] used the No-binding Theorem 
of Thomas-Fermi theory (see Lieb and Simon [2H|) instead of Baxter's estimate. 

If the fermions are described by the relativistic kinetic energy T Rcl or the magnetic 
Pauli operator T Pauh the proof of stability is not quite as simple. Note in particular that a 
potential that has a Coulomb singularity is not integrable to the fourth power and one can 
therefore not directly use the inequalities JSJ or The reader is referred to [35] and |2H] 
for the detailed proofs of Stability in these cases. 

5.1. MANY-BODY LOCALIZATION TECHNIQUES 

In many aspects of many-body theory the method of localizing has been a very useful 
technique. Dyson and Lenard used it heavily in their papers |12l I18| . 

We here discuss the two problems connected with localizing a charged system, i.e., the 
localization of the interaction and the localization of the kinetic energy. 

The simplest way to approach the localization of the kinetic energy is to do a Neumann 
localization. In fact, let — A^ u be the Neumann Laplacian for a unit cube centered at 
y &M 3 . Then in terms of quadratic forms defined on smooth compactly supported functions 
we may write 



where —A is the Laplacian on all of M 3 . Of course, there is a rescaled version to cubes of 
other sizes. The structure of this identity, i.e., that we have written the original operator as 
an integral over localized operators is characteristic for the way we shall write localizations 
in this section. The idea is then that for each value of the integration parameter one has a 
localized problem that one will estimate. Afterward the estimate is integrated. 

Another approach is to use a smooth localization. Let X be a C°° function with compact 
support. Let X y {x) — x(x — y). Then we have the identity 



which is often referred to as the IMS localization formula (see [B]). Instead of writing the 
localization as an integral one often uses a sum instead, but then one must introduce a 
partition of unity. 

For the problem of the charged Bose gas it turns out that we cannot use either of these 
localization methods. In this case we cannot completely localize the kinetic energy operator. 
We must still have some kinetic energy to control the variation between localized regions. 
The solution is to write the kinetic energy as a sum of a a high energy part and a low energy 
part. We then localize only the high energy part and use the low energy part to control the 
large distance variation. The result can be found in [231 . We shall state a simplified version 
here. 

Assume that X with < X < 1 is smooth has support in a unit cube centered at the 
origin in R 3 and that y/l — X 2 is also smooth. Let as before X y {x) = x(x — y). Let a*{y) 
be the creation operator for a particle in a normalized constant function in the unit cube 
centered at y € R 3 . Let V y be the orthogonal projection operator projecting on functions 




(10) 




(11) 
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orthogonal to constants in the cube. Then for any bounded set fl C M 3 and any < s < 1 
we have the following many-body kinetic energy localization estimate on the symmetric 
tensor product space L 2 (M 3 ) (the fully symmetric subspace of the full tensor product 

space) 

N „ N 



m 

3 



0) 
y 




+ J2 ( V a S(y + e > (2/ + e,) + l/2- ^aS(y)a (y) + 1/2) \dy 

-3vol(Q), (12) 

where ei, e2, e% is the standard basis in R 3 and e(x, s) — ► as s — > 0. Note that this is an 
explicit many-body bound in the sense that the right side of the operator inequality is not 
a one-body operator. The above estimate generalizes immediately to the situation where a 
charge variable is present. In this case there should be two creation operators, one for charge 
+1 and one for charge —1. There are rescaled version of (I12f) . which hold if we replace the 
unit cube by cubes of other sizes. 

The first sum on the right side of l|12l) is the localization of the high energy part of the 
kinetic energy. High energy means much larger than s~ 2 . 

The second sum on the right side of l|12l) represents the kinetic energy due to the vari- 
ations between localized regions. Using the Cauchy-Schwarz inequality one sees that for a 
state * e <g>^ m £ 2 (R 3 ) we have 

(*, {j^(y + e J )a a ( y + e ] ) + l/2- 

> (V K(y + ej)%(y + ej) + 1/2)*) - ^ (*, K(yK(y) + 1/2)*) 

This allows us to think of the last term on the right side of l|12|) as a discrete Laplacian 
acting on the function 

y ' * y/(*,(a* (y)a (y) + l/2)V). (13) 

For the Bose gas we shall conclude that most particles are in the constant function state 
(the zero momentum state). More precisely, this means that (^P, aj(y)a (?/)^') is almost the 
expected total number of particles in the unit cube centered at y. The expectation of the 
last term on the right side of Ijl2(l will therefore essentially give a contribution equal to the 
Laplacian of the square root of the density (assuming that we can approximate the discrete 
Laplacian by the continuous Laplacian). 

We finally come to the discussion of the localization of the interaction. For the Coulomb 
or Yukawa interaction this can be done using a method of Conlon, Lieb, and Yau [7]. Let 
as before X be a smooth function supported in the unit cube centered at the origin. There 
is an lo > depending on X such that for all fi > we have 

/x 2 ) ZiZjY^Xi - Xj) > / ^ X v {xi)Y^ +u (xi - Xj)x y (xj)dy - Nui. (14) 

^ ' l<i<j<N ^ R3 l<i<j<N 
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Note that the effect of the localization function x y on the right side is that for each fixed 
value of the integration parameter y only particles that live in the unit cube centered at y 
interact with each other. Again it is easy to see how this estimate changes under rescaling. 
A very elegant version of this estimate was given by Graf and Schencker |21| . 

6. THE LOWER BOUND IN DYSON'S CONJECTURE 

We shall now describe the main steps in the proof of the lower bound in Theorem 13. 21 The 
reader should look in |3()j for details. In each step we shall ignore certain errors and we shall 
not explain in details how these errors are estimated. In the detailed proof the errors are, 
in fact, only estimated at the very end when the main contributions have been identified. 

It is first of all important to understand that there are two relevant length scales in the 
problem. A long scale Lo, which is the diameter of the Bose gas and a short scale £q which 
is the distance on which the Bogolubov pairs interact. It will turn out that Lq ~ A -1 / 5 
and £ ~ A" 2 / 5 . 

Step. 1: We first localize the whole system into a large cube of size L ^> Lq ~ A _1 / D . 
On the other hand L should not be too large in order to allow us to control the volume 
error in l|12(l . This first localization is done using the Conlon-Lieb-Yau estimate (|14[) and an 
IMS-localization of the kinetic energy. Note that after rescaling the error in l|14(l will be 
Nuj/L < wA 6 / 5 . The IMS localization error (i.e., A / (Vx) 2 / / X 2 ) will be NL~ 2 < A~ 7 / 5 . 

Step. 2: We do a second localization into boxes of size £, where £q <C £ <C Lq. This 
localization is done using the many-body kinetic energy localization l|12f) together with the 
Conlon-Lieb-Yau estimate (|14f) . This time the error in using l|14|l is Ncu/£ <C Nuj/£q ~ 



The result of this localization is that the total ground state energy is estimated below 
by a sum of two terms. One term is the discrete Laplacian term discussed in the previous 
section. The other term is the local energy, which will eventually lead to the effective 
Hamiltonian for which we may apply the Bogolubov approximation. 

Step 3: Before proceeding with the analysis of the local short scale energy we introduce 
long and short distance cut-offs in the interaction. This is done as explained in the beginning 
of Section |5] 

Step 4: This is the final step in the reduction to the effective Hamiltonian. The localized 
two body potential is of the form W(x,x') — X y (x)V(x — x')x y (x'), where V is a cut-off 
Coulomb potential, i.e, V = Y tl — Y v for some appropriate /i -C v. Here X y has support in 
a cube of side £ centered at y. 

Let a* az be the creation operator for a particle of charge z = ±1 in a state u a , a = 0, 1, . . ., 
where {u a } is an orthonormal basis for the Hilbert space of L 2 on the cube of size £ centered 
at y and uq is the constant function. It is not otherwise important how this basis is chosen. 
We may write the full two body interaction as 



The fourth step in the proof is to show that one may ignore (for a lower bound) all terms in 
this sum which do not have precisely two of the parameters a, P, 7, S equal to zero. Moreover 
one may also ignore terms for which a = 7 = 0or/3 = <5 = 0. This step in the proof is rather 
technical and uses heavily that we have been able to introduce cut-offs in the potential. 



luN 7 / 5 . 




z.z' — ^zl a/376 
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The resulting two body interaction may be written 

/ ^3 / V(p) ^ ' ZZ '(Kp,z b +p.z> + b *-p.z b -p,z> + Kp.z b -p.z> + b +p,z b -p,z') d P' 
2{27T) J P£ R3 z,z'=±l 

where V denotes the Fourier transform of the potential V, denote the total number of 
positively and negatively charged particles in the cube respectively, and 

K,z = (e 3 vT 1/2 <CP y X y e i naoz, 

where a*(/) creates a particle of charge z in the state / and as before V y is the projection 
orthogonal to constants. It is easy to see that for fixed p the operators b± p z satisfy the 
conditions in Theorem ED 011 the Foldy-Bogolubov method. 

It is also easy to see that the localized kinetic energy resulting from l|12|) may be estimated 
below by 

7777^ f *(*>) K P .z b T P ,zdp, 

where t(p) = ^£ 3 p i {p 2 + £s~ 2 )-\ 

It is now immediate from the Foldy-Bogolubov method in Theorem 14.31 that the local 
ground state energy is estimated below by 



"97^3 / (*(p) + "v^(P)) ^ \/(t(p) 2 + nP>yV(p)dp 

where v y = v + + v~ (which depends on the location of the cube, although this had been 
suppressed for v^). If we now ignore the various cut-offs (which may be justified) and simply 
replace t(p) by \£ z p 2 and V{p) by A-k /p 2 we see that the above integral gives — Jz^ 4 i?~ 3 / 4 = 
-J(u y /e 3 ) 5 / 4 e 3 , with J defined in Theorem EH1 

Step 5: In order to control the errors encountered in reducing to the effective Hamilto- 
nian as well as in treating the discrete Laplacian term we must conclude that most particles 
are in the condensate, as explained in the previous section. To do this, we shall localize the 
number of excited particles. To be more precise, we may think of the Hilbert space as being 
a direct sum over subspaces with a definite number of particles in the condensate. The 
Hamiltonian is not diagonal in this representation. We want to conclude however that re- 
stricting any given state to only a small finite number of the subspaces will not significantly 
change its expected energy. This is achieved using the following result from |29| . 

Theorem 6.1. (Localization of large matrices) 

Suppose that A is an N x N Hermitian matrix and let A k , with k — 0, 1, . . . , N — 1, denote 
the matrix consisting of the k th supra- and infra- diagonal of A. Let ip S C N be a normalized 
vector and set dk = (i/j,A k ip) and A = (ip,Aip) = dk- (ip need not be an eigenvector 

of A.) 

Choose some positive integer M < N . Then, with M fixed, there is some n £ [0, N — M] 
and some normalized vector <f> £ C N with the property that <j>j = unless n+l<j<n + M 
(i.e., cj) has length M) and such that 

„ M-l JV-l 

W^<A + ^5> 2 K| + c]L |d fc | , ( 15 ) 

k=l k=M 
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where C > is a universal constant. ( Note that the first sum starts with k = 1 .) 

Step 6: The final step is to combine the two parts of the energy described in Step. 2. The 
local energy was (approximately) bounded below by — J(t , y /^ 3 ) 5/ ' 4 £ 3 which when integrated 
over y and normalized by the volume of the cube gives 

-J [ (v y /l 3 f 4 dy. 
Jr 3 

The other part of the energy is essentially the kinetic energy of the function in (where 
we had actually ignored the charged variable). Using the result of Step 5 that most particles 
are in the condensate we may write this kinetic energy as (approximately) 




If we use that the total number of particles is N we have the condition that L 3 v y /l 3 dy = N. 
Thus if we define <j>(y) — N~ 4 ^ 5 v N -i/n y /£ 3 we see by a straightforward scaling argument 
that J 4>{x) 2 dx = 1 and that the ground state energy is approximately bounded below by 

N 7/5 (lf (V4>(x)fdx-J f {4>(x)f 2 dy\ 
which is of course bounded below by minimizing over all as in |7|. 



7. THE UPPER BOUND IN DYSON'S CONJECTURE 

To prove an upper bound on of the form given in Dyson's conjecture Theorem 13.21 we 
shall construct a trial function using as an input a minimizer (f> for the variational problem 
on the right side of (Q. That minimizers exist can be seen using spherical decreasing 
rearrangements. Define 4>q(x) — A r3 / 10 (/)(7V 1 / 5 a;). 

Let cf> a , a = 1, ... be an orthonormal family of real functions all orthogonal to <pQ. We 
choose these functions below. 

We follow Dyson and first choose a trial function which does not have a specified 
particle number, i.e., a state in the space ®^ =0 Hn, {7~(-o = C) more precisely in the bosonic 
subspace, i.e., the bosonic Fock space. We shall evaluate the expected value of the operator 
Sat=o in our state. 

As our trial many-body wave function we now choose 

*= H^-Wex* (-Ao+AooS + + AooS_- ]T £ ^'<X,*' ) |0) , (16) 

where a* z is the creation of a particle of charge z = ±1 in the state <p a , |0) is the vacuum 
state, and the coefficients Ao, Ai, . . . will be chosen below satisfying < A Q < 1 for a ^ 0. 
It is straightforward to check that \E' is a normalized function. 

Dyson used a very similar trial state in |llj . but in his case the exponent was a purely 
quadratic expression in creation operators, whereas the one used here is only quadratic in 
the creation operators a* z , with a ^ and linear in Oq_|_. As a consequence our state will 
be more sharply localized around the mean of the particle number. 



February 7, 2008 14:39 WSPC/Trim Size: lOin x 7in for Proceedings solovej-icmp-june3 



15 



Consider the operator 



°° \ 2 



(17) 



a=l 

A straightforward calculation of the energy expectation in the state \1/ gives that 

( *, jr H N *\ = A 2 , /"(V0o) 2 + Tr (Tr) + 2A 2 Tr (iC (v - y/T(T+l) 
\ w=o / J 

where K is the operator with integral kernel 

y) = (f> (x)\x - y\~ Vo(y)- 

Moreover, the expected particle number in the state \E r is 2A 2 , + Tr(r). In order for \& to be 
well defined by the formula (jTI))) we must require this expectation to be finite. 

Instead of making explicit choices for the individual functions <p a and the coefficients 
A Q , a/0we may equivalently choose the operator T. In defining T we use the method of 
coherent states. Let X be a non-negative real and smooth function supported in the unit ball 
in M 3 , with J X 2 = 1. Let as before N~ 2 / 5 < £ < N' 1 ^ and define Xt{x) = £~ 3 / 2 x(x/ '£). 
We choose 

T(x,y) = (2tt)- 3 f f(u,\p\)V£ \6 u , p ){6 u , p \Vldudp 

where V^ Q is the projection orthogonal to <fio, 

Ou, P (x) = exp(ipx)Xi(x - u), 

and 

p 4 + 167rA^ (u) 2 \ 



/(«, N) 



2 V(P 4 + 32^A20 O ( U )2) 1 / 2 



We note that T is a positive trace class operator, T(j>o = 0, and that all eigenfunctions of 
r may be chosen real. These are precisely the requirements needed in order for L to define 
the orthonormal family <fi a and the coefficients A a for a ^ 0. 

We use the following version of the Berezin-Lieb inequality jJJES]- Assume that is 
an operator concave function on M. + U{0} with £(0) > 0. Then if Y is a positive semi-definite 
operator we have 

Tr (Y"£(L)) > (2^)- 3 f e(/(«, \p\)) (e u , P ,V£ YV£ o e UtP ) dudp. (18) 



We use this for the function £(f) = y/t(t + 1). If Y = I then ljT%|) holds for all concave 
function £ with £(0) > 0. Of course, if £ is the identity function then H18|l is an identity. 

This reduces proving an upper bound on the energy expectation (|17fl to the calculations 
of explicit integrals. After estimating these integrals one arrives at the leading contribution 
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(in N) 

Ao J Wo) 2 

+ JJ Q/ + 2A^oH 2 ^) /(«, M) - ^2A^o(«) V/(«, bl) + 1) dpd« 

where J is as in (JJJ. 

If we choose Ao = y/N/2 we get after a simple rescaling that the energy above is N 7 ^ 5 
times the right side of Q (recall that <f> was chosen as the minimizer). We also note that 
the expected number of particles is 

2A 2 , + Tr(r) = N + 0(N 3/5 ), 

as N — > oo. 

The only remaining problem is to show how a similar energy could be achieved with a 
wave function with a fixed number of particles N, i.e., how to show that we really have an 
upper bound on En- We indicate this fairly simple argument here. 

We construct a trial function Vf 1 ' as above, but with an expected particle number N' 
chosen appropriately close to but slightly smaller than N. Using that we have a good 
lower bound on the energy En for all N we may, without changing the energy expectation 
significantly, replace \E ,/ by a normalized wave function \1/ that only has particle numbers 
less than N. Since the function N i— > Ejy is a decreasing function we see that the energy 
expectation in the state \& is, in fact, an upper bound to En- 
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